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Understanding the stochastic dynamics of tracer particles in active fluids is important for iden- 
tifying the physical properties of flow generating objects such as colloids, bacteria or algae. Here, 
we study both analytically and numerically the scattering of a tracer particle in different types of 
time-dependent, hydrodynainic flow fields. Specifically, we compare the tracer motion induced by 
an externally driven colloid with the one generated by various self-motile, multi-sphere swimmers. 
Our results suggest that force-free swimmers generically induce loop-shaped tracer trajectories. The 
specific topological structure of these loops is determined by the hydrodynamic properties of the 

J.—S , microswimmer. Quantitative estimates for typical experimental conditions imply that the loops 

»vj ' survive on average even if Brownian motion effects are taken into account. 

^ PACS numbers: 47.63.Gd, 05.40.-a, 47.63.mf 

(N 

,"! I. INTRODUCTION 

o 

^ ' Scattering processes are ubiquitous in Nature, ranging from elementary particle collisions at subatomic scales to 

gravitational encounters in galactic clusters. Traditionally, scattering experiments have played an important role in 
elucidating the interactions between 'non-living' physical objects. Nowadays, modern experimental techniques allow 
us to track the motion of individual microorganisms HHj] ^as beautifully illustrated by recent high-speed microscopy 
observations of Volvox [5] and Chlamydomonas reinhardtii [6|,|7|. These and similar experiments on colloidal systems [8| 

^ ' suggest that it should be possible in the near future to systematically study the effective hydrodynamic interaction 

O ! forces generated by algae, bacteria or artificial microswimmers |9l.ll0l| through suitably designed biophysical scattering 

O^ ' experiments. 

Additional motivation for studying swimmer-tracer scattering comes from recent experiments by Leptos et al. [y|, 

^ . who observed that a passive tracer particle exhibits non-Gaussian diffusion when surrounded by a dilute suspension 

^ ' of sclf-motilc Chlamydomonas reinhardtii algae. Qualitatively, the anomalous diffusive behavior arises because the 
^D , algae modify the velocity field of the fluid, thereby occasionally accelerating the tracer particle to relatively large 
^^ ■ velocities. A satisfactory quantitative explanation for this phenomenon is still missing, and a better understanding 
of the underlying elementary scattering processes is a key step en route to deriving the corresponding generalized 
diffusion equation. 
C^ ' The present paper aims to identify generic features that should be observable when a small colloidal tracer particle 

is scattered by an artificial or natural microswimmer. Specifically, we address the following problems: How does the 
tracer dynamics differ depending on whether a swimmer's motion is externally forced or internally generated? Is it 
sufficient to consider an effective time-averaged description in order to correctly predict the tracer trajectories - or 
is it necessary to use a time-resolved model that captures the details of the swimming stroke? How does Brownian 
motion, which plays a non-neglible role for (sub-)micron-sized colloidal particles, affect the tracer dynamics? To 

; ^ ' answer these questions, we will compare both analytically and numerically the dynamics of a tracer in the presence 

Ci ' of various simplified model swimmers. 

The first part of our discussion focusses on deterministic scattering processes. In order to clarify how the effective 
interaction range determines the asymptotic state of the tracer, we shall consider a simplified power-law interaction 
model and, as a more realistic example, tracer displacement due to an externally driven colloid. The main part of 
the paper is dedicated to the scattering of a tracer by 2-sphere and 3-sphere swimmers that generate predominantly 
dipolar and quadrupolar fiow fields, respectively. We show that, asymptotically, the tracer motion in the presence of 
a self-motile, force-free swimmer converges to a closed loop. The shape and direction of the loops is a signature of 
the specific properties of a given swimmer. For experimentally accessible parameters, Brownian motion effects may 
dominate the trajectories of an individual tracer particle. In the last part of the paper, we will demonstrate that, after 
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averaging over a few hundred to thousand samples, the mean tracer trajectories look very similar to those obtained 
in the deterministic limit. 

II. MATHEMATICAL MODEL 

We consider the motion x{t) of a passive, colloidal tracer particle in a fluid due to the presence of an active object 
('swimmer'). The latter, which may be an externally forced colloid, an artificial microswimmer or a micororganism 
(e.g., an alga or a bacterium), is described by a phase space coordinate vector T{t). For example, for a spherical 
colloid with position vector X{t) and velocity X{t) := dX(t)/dt, T{t) is simply given by T{t) — {X{t),X{t)). For 
more complex objects, such as multi-sphere swimmers jlll - ll3l |. F comprises all the coordinates necessary to uniquely 
specify the motion of the swimmer. Throughout, we shall assume that, in good approximation, the tracer particle 
does not affect the swimmer motion, i.e., F(<) is taken to be independent of x{t). 

A. Langevin dynamics of the tracer 

In principle, when studying the scattering of a tracer particle by a self-swimming object, one can distinguish (at 
least) two different approaches: 

(i) analyzing the approximate mean motion of the tracer particle in the time-averaged (or stroke-averaged) effective 
flow field of the swimmer. 

(ii) computing the exact, time-resolved motion of the tracer particle in the full oscillatory flow field of the swimmer. 

Below we will compare results from both methods for multi-sphere swimmer models [121 |l3| . We will start from the 
time-averaged description (i) which is less accurate, but allows us to obtain analytical estimates. When adopting this 
coarse-grained approach, the swimmer state is approximated by F(t) = {X(t), X{t)), and the swimmer is assumed to 
move at constant velocity X{t) = V along the straight line 

X{t)^Xo + tV. (1) 

The effective stroke-averaged swimmer velocity V results from the microscopic swimmer dynamics. 

Let us assume we know the velocity field v{x\T{t)) that is generated by a given swimmer configuration F(t). In 
this case, we can model the motion of a micron-sized tracer particle by the overdamped Langevin equation 

±{t) - v{x{t)\r{t)) + {2Do)^'H{t). (2a) 

This equation is valid in the Stokes (zero Reynolds number) regime, which is realized to good approximation under 
typical experimental conditions [a-0|- The second contribution on the rhs of Eq. ((2a)) is Gaussian white noise £,{t) = 
(^i(t)), which describes thermal fluctuations in the fluid and is characterized by 

(6(t))=0, {^.mj{s))=5.,,5{t-s). (2b) 

For a spherical tracer particle of radius oq, the thermal diffusion constant Dq is given by the Stokes formula 

Do = kTho = fcT/(67rr;ao), (2c) 

where oq is the tracer radius and 70 the Stokes friction coefficient. For example, considering oq — 1/xm and water 
at room temperature T = 25°C with viscosity rj — lO^^kgm^^s^^ and kT — 4.11 x 10~^^kgm^s~^, we have Dq = 
0.22^m2s"i. 

B. Initial conditions 

We intend to analyze swimmer-tracer scattering processes for well-defined, reproducible initial conditions. For this 
purpose, we choose the coordinate system such that the spherical tracer particle is initially located at the origin, 
x{0) — 0. The swimmer (characteristic length scale A) starts at Xq = {Xo,Yo,0) with lYol > ao + A and moves in 
the (z = 0)-plane with an average velocity V = (V, 0, 0) parallel to the x-axis, see Fig. [TJ 
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FIG. 1: Swimmer-tracer scattering. The swimmer (S) starts at Xq — {Xo,Yo,0) and moves with velocity V = (V^O, 0) parallel 
to the a;-axis. The trajectory of the tracer (T), initially located at a;(0) — (0,0,0), depends on the flow field generated by the 
swimmer. 

III. TRACER SCATTERING IN TIME-DEPENDENT FLOW FIELDS 

Our main objective is to understand the motion of the tracer in the velocity field v{x\T(t)) of an externally driven 
or self-motile object. In order to simplify the notation, it is convenient to define the distance vector between tracer 
and swimmer by 

r{t)^x{t)~X{t), (3a) 

and to introduce two unit vectors 

with n{t) characterizing the orientation of the swimmer motion. In the examples considered below, the effective 
hydrodynamic interaction between swimmer and tracer is mediated by non-diagonal tensors iif — H{\r\,r,h) = (Hij), 
which relate the fluid velocity at the position x{t) to the velocity V of the swimmer through 

V = HV. (4) 

The exact structure of H encodes the details of the propulsion mechanism. For example, consider an externally 
forced colloid of radius A that is pulled at constant velocity V. If the colloid passes the tracer particle at a sufficiently 
large distance \r{t)\ := \x{t) — X{t)\ ^ A + qq, then H is given by the Oseen tensor 

3A 

Hij ^ ^AStj + hf-j) . (5) 

By contrast, the effective velocity field generated by a self-swimming microorganism usually decays more rapidly 
as l^l"", a > 2. In Sections IIII Bl IIII CI and IIII Dl we shall compare the scattering of tracers by externally forced 
colloids and multi-sphere swimmers. However, in order to demonstrate how the effective interaction range a affects 
the scattering process, it is instructive to analyze a simplified power law interaction model first. 

A. Simplified power law interaction (a-model) 

Let us assume that a hypothetical 'swimmer' of radius A generates a fluid velocity field of the simple power law 
form 

v{x\Xit),V) = V. (^)". (6) 

The dimensionless parameter k quantifies the interaction strength. In order to understand the basic dynamics of a 
tracer in this field, we will focus on the deterministic limit Dq = in the Langevin equation (Pa|l . Brownian motion 
effects (corresponding Dq > 0) will be discussed later for more realistic flow fields. 

Adopting the initial conditions V — {V, 0, 0), X{0) — {Xq, Yq, 0), and a;(0) = as depicted in Fig. [TJ the swimmer 
moves in the x-direction and its x-position at time t is given by X{t) — Xq + tV . The isotropic power law model ^ 



corresponds to the diagonal tensor Hij = K{A/\r\)°'Sij, which imphes that the tracer particle also moves along the 
X-axis, i.e., x{t) = (x(i),0, 0) where 



x{t) ^ Vk 



{[X{t)-x{tW+Y,^V/^ 



(7) 



It is possible to find exact, implicit solutions of Eq. ([7]) for a — 1,2, .. . by transforming to the comoving swimmer 
frame. For example, for long-range interactions with a = 1, we obtain, in the original coordinate system. 
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The exact results (|S]) are transcendental equations for a;(i) which can be numerically inverted for a given value t > 
to give the trajectory of the tracer particle. In particular, these equations cover both the near- and far-field scattering 
behavior. However, to gain additional analytic insights, it is useful to consider in more detail the far-field scattering, 
corresponding to \x(t)\ ^ max{|yo|, I-'^IOI} ^'^'^ k'^~"A^ ^ Y^ . In this limit one obtains, for a = 1, the approximate 
result 



x(i) ~ kA\o% 



x{t) + ^xWTy^ 



Inserting X{t) = Xq + tV, we see that asymptotically (for t — )■ oo) 

X[t) ~ KAlOg 



-'^O + \/Xq + ^0 



(9a) 



(9b) 



Thus, for V ex |r| ^ the solution is logarithmically divergent in time, implying that the tracer particle is slowly dragged 
away to infinity. Defining the asymptotic displacement as a function of a by 
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(10) 



the divergence implies (Ax) 5" = oo. Performing a similar analysis for a = 2 yields 

A 
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\Yn 



arctan I ^V -tan fl^ 



and thus, in the limit t — >■ oo. 
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i.e., for a = 2 the tracer particle comes to rest after a finite distance (Ax)^ . Analogous results can be obtained for 
arbitrary values a > 2. For example, for a = 3 
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FIG. 2: Deterministic tracer trajectories for three different model interactions, (a) Tracer displacement in the direction of 
the swimmer's motion as a function of time t for the simplified a-model with k = 1 and the Oseen interaction, (b) Tracer 
trajectories in the {x,y)-p\a.ne. For colloidal-type long-range interactions that scale with distance as \r\~^ the tracer is dragged 
to infinity - the divergence is logarithmic in time. By contrast, for shorter range interactions with |rj~",Q: > 2 the maximal 
tracer displacement remains finite, see red curve/triangles. Solid lines correspond to the exact result ^, and dashed lines 
to the approximate far-field result ((5)). The symbols were obtained by numerically integrating the deterministic equations of 
motion (|2aP with Do — using a simulation time step At = 2.2 x 10~^s. Parameters: swimmer velocity V — 100^/s, swimmer 
radius A = 10/im, tracer radius a — l/xm, impact parameter Yo = 10^, initial swimmer position Xo — —25^4. 



To briefly summarize, for long-range hydrodynamic interactions with v oc \r\~^ the tracer experiences an infinite 
displacement, whereas for shorter range interactions with v ex |r|~",a > 2 the displacement is finite. These results 
are illustrated graphically in Fig. [H 

The simple power law model (jG)) is useful for illustrating how the range of the effective hydrodynamic interactions 
affects the asymptotic state of the tracer particle. However, this truncated model neglects fiuid velocity components 
that act perpendicular to the swimmer's direction of motion. In the remainder, we will focus on more physical models 
that account for these effects. As the first example, we consider the scattering of a tracer by a driven colloid. 



B. Externally forced colloid 



Consider a spherical colloid (radius A) dragged at constant velocity V past a smaller tracer particle (e.g., by a 
movable laser trap). If the distance between tracer and colloid is sufficiently large, \r{t)\ :^ A + qq, the velocity field 
components Vi = HijVj experienced by the tracer are determined by the Oseen tensor ([5]) and we can approximate 
r{t) ~ X{t). Considering the deterministic limit Dq = and the same initial conditions as before, we then find that 
the trajectory of the tracer particle in the (z = 0)-plane is determined by the differential equation 
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Owing to the non-vanishing off-diagonal components of the Oseen tensor, the tracer experiences a transverse force 
(velocity field) component in the y-direction. Equation (|13p can be directly integrated, yielding 



x{t) = kA{ 2 log 
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where X{t) = Xq + tV . Letting i — > oo, we find that asymptotically 
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FIG. 3; Simple multi-sphere swim.m.er m,odels. (a) The 2-sphere swimmer with variable radius spheres generates a dipolar 
net flow field that can be either extensile or contractile with respect to the swimmer axis, see Figs. U (b) The 3-sphere 
swimmer [I2| generates a flow fleld that is essentially quadrupolar [Q, see Figs. [T] Consequently, the trajectories of tracer 
particles are qualitatively different when scattered by the two swimmers. 

Note that the logarithmic divergence predicted by Eq. (|15a|) is in agreement with the earlier result (jllap . i.e. the 
tracer particle slowly follows the colloid to infinity. 

The deterministic result p^ can be expected to correctly describe a single tracer trajectory only if Brownian fluc- 
tuations are negligible corresponding to the low-temperature, high-viscosity regime. For a micron-sized tracer particle 
in water at room temperature, Brownian motion effects can become relevant if the speed V of the 'swimmer' colloid is 
too low and/or the distance between tracer and swimmer colloid becomes too large. In this case, Eq. (J14p provides an 
approximate description of the ensemble mean when the scattering experiment is repeated using identical (i.e., deter- 
ministic) initial conditions. The statistical fluctuations around the mean trajectory P^ can be estimated by ^/2Dot, 
where Dq ~ 0.22/^m^/s for a micron-sized tracer. This suggests that for the parameters and initial conditions used in 
Fig. [5J the characteristic features of the solution p^ should remain observable even for a single tracer trajectory. 

Having discussed the scattering of a tracer by an externally driven colloid, we shall focus on self-motile swimmers 
in the remainder of the paper J12h17| |. In particular, we would like to identify characteristic features of the tracer 
trajectories that provide a signature of truly self-swimming objects, which do not generate a net force during a 
swimming stroke. Specifically, we are interested in understanding how the tracer trajectories differ depending on 

(i) whether a force-free swimmer is extensile ('pusher') or contractile ('puller') ; 

(ii) whether a force-free swimmer generates a dipolar or a quadrupolar flow field. 

For this purpose, we will study simple multi-sphere swimmers that are analytically tractable and can be easily tuned 
from extensile to contractile. An additional advantage of the simplified models is that the und erly ing microscopic 
swimming mechanisms can be implemented numerically using the methods described in Ref. J16l Il7| . This allows us 
to compare analytic estimates, which are based on a time-averaged description, with numerical results for the exact 
time-resolved dynamics. 

C. Dipolar two-sphere swimmer 

As the first example of tracer scattering by a self-motile object, we consider a slender 2-sphere swimmer consisting 
of two beads (radii A1.2) at positions Xi,2{t) that are connected by an oscillating leg of length 

L(t)=£ + ^sm{ujt + 9) , £>C>0, (16) 

see Fig. |3] (a). The orientation vector of the swimmer is defined by 



If the sphere radii Ai and A2 are constant then this 2-sphere swimmer reduces to an oscillating dumbbell ^^ [T^, [ij] , 
which is prevented from swimming by the scallop theorem |15l. Il8l|. A simple, self-motile 2-sphere swimmer can be 
obtained if the size of the spheres changes periodically in time according to 

Ai{t) =A + Xsm{ujt + ,pi) , A2it) =A + Xsm{ujt + (j>2) , A > A > 0. (18) 

For simplicity, we neglect secondary flow contributions produced by the expansion and shrinkage of the spheres by 
assuming that (to good approximation) they merely change their hydrodynamic resistance in a manner consistent 
with a changed radius. Considering the limit A,^ <^ A,i and averaging over a stroke-period [t,t + 27r/a;], one can 
show that the 2-swimmer moves at an average velocity V — Vh, where 
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^ - 4^A ( ^ 1 ( 1 - — ) sin(A</.) cos(A0) (19a) 



and 



A(/):=^(0i-02), A0:=i(0i+02)-e. (19b) 

The derivation of Eq. (fT9|) is analogous to those for other multi-sphere swimmers [ij, [l^ . By choosing the phase 
parameters (f>i,(j)2, and 9 such that A(j> E (0,7r) and A9 e (— 7r/2,7r/2), we have V > 0, i.e., the swimmer moves in 
the direction of positive n. 

The stroke-averaged velocity field generated by the swimmer at large distances is obtained by summing the Stokeslets 
of the two spheres, performing a far-field Taylor (multipole) expansion and averaging the leading far-field contributions 
over a swimming stroke [13l] . By means of this procedure, one finds that the averaged far- field fluid flow generated by 
the swimmer is dipolar, 

v{x{t)\X{t),V) = Vk2(^j^^ [3(nf)2-l]f, (20a) 

with unit vectors n and f as defined in Eq. (j3bp and 
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4:\AJ \ 2£ J tan(A</.) 

Note that, depending on the choice of Acf) and A9, the coefficient K2 can be either negative or positive for a swimmer 
with V > 0. Swimmers with K2 > are extensile 'pushers', whereas negative values K2 < correspond to contractile 
'pullers'. In our simulations, we used 6* = 0, Acf) — A9 — 7r/4 to realize an extensile swimmer and 9 = 0, Acf) = 
37r/4, A9 = 7r/4 in the contractile case. The structure of the dipolar flow fields (I^Ul) for these parameter values is 
illustrated in Figs. 21 The arrows show the normalized velocity field i) := i'/|'u|- The color shading indicates the 
projection (ri)), i.e., dark (bright) areas correspond to regions where the mean fiow points towards (away from) the 
swimmer. 

Tracer motion in the deterministic limit.- Based on Eqs. (1201) . it is possible to obtain an analytic estimate for the 
mean tracer trajectory in the deterministic, far-field limit. To this end, we consider initial conditions with n = (1, 0, 0) 
as depicted in Fig. [T] and approximate 

r{t)^X{t), f,v~ ^^* 1 (21) 

,/xWTYi' ^ ^ 

with X{t) = Xq + tV denoting the swimmer's x-coordinate. The approximations (PT|) hold in the far-field scattering 
limit where Eq. pO)) is valid. Integrating the equations of motion ((2a|) for Do = yields the trajectory x{t) = 
{x{t),y{t), 0) of the tracer particle in the (z = 0)-plane 

^m ~ _„,A!4(Hi±2iL)_iHW!±22\ ,22.) 
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FIG. 4; Dipolar fluid flows generated by 2-sphere swimmer.- The diagrams show the stroke-averaged fluid velocity field (|20p 
for two different swimmer setups: (a) extensile pusher with K2 > and (b) contractile puller with K2 < 0. In both cases 
the swimmer is located at the origin, points along the a;-axis, and swims in the positive x-direction. Swimmer parameters: 
£ = 20AJm, ^ = A = 0.1^, A = 0.1A,lj = 5 x lO^s"^ with A(j) ^ AO = ir/4 for extensile swimmers and A(j) = 37r/4, AO = 7r/4 
in the contractile case, yielding V ~ 14.7/im/s and H2 ~ ±8.8. The flow flelds are rotationally invariant about the x-axis. 
Arrows indicate the local direction v := v/\v\ of the stroke-averaged flow (|20p . Dark (bright) areas show the projection (rv), 
corresponding to regions where the mean flow points towards (away from) the swimmer. 



Letting t — > CXI, we see that asymptotically 
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According to Eqs. (P^. the tracer particle is displaced by a finite distance reflecting the fact that the effective interac- 
tion is short-range a > 2. Furthermore, Eq. (j23l) predicts that, if the swimmer starts very far away, corresponding to 
the limit Xq — >■ — c», the tracer returns to its initial condition a;(0) = 0, thus forming a closed loop in the (z = 0)-plane. 

This is illustrated in Fig. [51 which shows tracer trajectories in the {x, j/)-plane. The dotted lines in both diagrams 
indicate the analytic estimates from Eq. (|22l) . Symbols in Fig.[S](a) were obtained by numerically integrating Eq. (Pa|) 
for Dq — for the time-averaged dipolar velocity field ([20)) using the same swimmer parameters as in Fig. S] For 
comparison, we also plot in Fig.[5](b) the result from simulations that resolve the microscopic dynamics of the swimmer 
(details of the numerical implementation are given in Sec. IIV[) . In the latter case, one finds that the tracer particle 
performs a small oscillator motion around its mean trajectory. Hence, by comparing Figs. [5] (a) and (b), one readily 
observes that the stroke-averaged description is able to capture the main features of the mean particle trajectory. 

Brownian motion effects.- As exemplified by Eq. (j20p . the effective hydrodynamic flow flelds of self-swimming 
objects decay more rapidly with distance than those of externally forced colloids. Hence, the hydrodynamic displace- 
ment of a tracer particle in the far-field fiow of a force-free swimmer will, in general, be smaller, so that thermal 
Brownian motion of the tracer becomes more relevant. This raises the question as to whether loop-like structures 
could actually be observed under typical experimental conditions, e.g., similar to those considered by Leptos et al. @ 
who tracked the motions of algae in water at room temperature. In order to estimate how Brownian effects may affect 
the loop-like structures, we numerically integrated the stochastic equations of motions (j2ap for both the time-averaged 
flow fleld ([20]) and the time-resolved swimmer dynamics, using a diffusion constant Dq — 0.22/im^/s corresponding to 
a tracer particle of radius oq = 1/im and water at room temperature. 

Figure [B] (a) depicts the mean tracer trajectories obtained numerically for the time-averaged model for an extensile 
(red triangles) and a contractile swimmer (blue circles). Figure [6] (b) shows the corresponding results for the time- 
resolved swimmer dynamics. The numerical data points (flUed symbols) in Figs. [S] were obtained by averaging over 
1000 sample trajectories with identical initial conditions. In both cases, we used the same swimmer parameters as in 
Figs.[5l but the swimmer was started closer to the tracer at Xq = — 35A. The simulation time per run was limited to 
10s. For comparison, isolated Chlamydomonas algae typically swim on an almost straight line for about 5 — 10s Q 
before changing their directions due to tumbling and rotational Brownian motion. The statistical error bars in Figs.[S] 
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FIG. 5: Deterministic tracer scattering by a dipolar 2-sphere swimmer. In the limit Do = the tracer trajectories converge 
to closed loops if the force-free 2-sphere swimmer starts sufficiently far away from the tracer (arrows indicate how the loop is 
traversed). The dotted lines in either diagram indicate the analytic estimate from Eq. ()22|) . (a) Results for the stroke-averaged 
dynamics. Symbols were obtained by numerically integrating Eq. H2ap for Do ~ for the dipolar flow field H20|) using the same 
parameters as in Fig. U (b) Results for the time-resolved dynamics. Symbols were obtained from numerical simulations of the 
tracer trajectories as described in Sec. IIVI The inset shows the details of the oscillatory motion of the tracer particle near the 
turning point at the top-right corner of the extensile (right hand) path. The plots are based on initial conditions Xo = — 350A, 
Yo = IQA for the swimmer and a;(0) = for the tracer. The tracer trajectories are depicted for the time interval [0, 100s]. The 
swimmer parameters are the same as those in Fig. [J] Tracer radius: ao = Ifim. Simulation time step: At = 2.2 x 10~^s. A 
comparison of the two diagrams implies that the stroke-averaged description is able to capture the main features of the mean 
particle trajectory. 
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FIG. 6: Brownian tracer scattering by a dipolar 2-sphere swimmer. For Do > loop-like scattering patterns can be recovered 
by averaging over many sample trajectories with identical initial conditions a;(0) — 0. (a) Results for the stroke-averaged 
dynamics obtained by numerically integrating Eq. (|2ap with Do — 0.22/im /s. (b) Results from the corresponding time- 
resolved simulations. The filled symbols in both figures were obtained by averaging over 1000 sample trajectories and dotted 
lines represent the analytic estimate from Eq. (|22|) . Statistical error bars correspond to the sample standard deviation, divided 
by 10 X \/i for better visibility. Prior to rescaling the length of the error bars agreed with the theoretically expected value \/2Dot. 
Compared to Fig. [5l we used the same swimmer parameters but the swimmer was started nearer to the tracer at Xq — —10^4 
and sample trajectories were recorded for the shorter interval t G [0, 10s]. This explains why the loops are not fully closed. 



correspond to the sample standard deviation and were divided by 10 x \/t for better visibility. Prior to this rescaling 
the length of the error bars agrees well with the theoretically expected value \/2Dat. 

To summarize, the diagrams in Fig. [5] suggest that after averaging over a sufficiently large sample size the loop-like 
patterns can be recovered. Individual trajectories may however look very different due to Brownian motion. 



D. (Quasi-)quadrupolar three-sphere swimmer 

As the final example, we consider the scattering of a tracer by a linear 3-sphere swimmer [12l Il3l|. as illustrated in 
Fig- El (b)- The main difference when compared to the 2-sphere swimmer occurs because the flow field of this swimmer 
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FIG. 7: Quasi- quadrupolar fluid flows generated by 3-sphere swimmers. Structure of the stroke-averaged flow field (|27p for three 
different parameter sets: (a) extensile pusher with 0.21 — ^r > £,f ~ O.li, (b) symmetric swimmer with ^r = ^f ~ \/0.02^, 
and (c) contractile puller with 0.1^ ~ £,r < £,f ~ 0.21. In all three cases, the swimmer is located at the origin, points along the 
a;-axis, and swims in the positive a;-direction. The flow fields are rotationally invariant about the a::-axis. Arrows indicate the 
local direction i) :— v/\v\ of the stroke-averaged flow field (|27[) . Dark (bright) areas show the projection (rv), corresponding to 
regions where the mean flow points towards (away from) the swimmer. Note that in all three cases the quadrupolar contribution 



dominates at intermediate distances. 
V ~ 12/im/s. 



Swimmer parameters: £ = 10fj,ra, A = 0.2£,uj = lO'^s ^,Acj> = it/2, corresponding to 



is essentially quadrupolar, see Fig. [T] The lengths of the rear leg (i?) and the front leg (F) of the 3-sphere swimmer 
vary in time as 



LF,R{t) = i + ^F.R s[n{ujt + (i>F,R) 



n > iR.F > 0, ^ > A > 0, 



(24) 



where t is the mean leg length and i,R,F denotes the oscillation amplitude. If the initial phases <i>R,F are chosen such 
that 



A(/) := (/>R 
then the swimmer moves with stroke-averaged speed 

7 



>0, 



F=-A^sin(A0)(^ 



(25) 



(26) 



in the direction of the front leg. The time averaged flow field at distances far away from the 3-sphere swimmer is [ij 



v{x[t)\X(t),V) = Vn2 



Vk^ 



A 



\r{t)\ 
A 



[i{hrf - l] r + 

{ [lb{hrf - 9(nf)] r - [2,{hrf - l] n] 



(27) 



with unit vectors n and r as defined in Eq. (j3b[) . The constants K2.3 are determined by the swimmer's geometry and 
read 



K.2 



9 a 



2 c2 



32 Ai 



51 
"^^56U 



(28) 



The quadrupole coefficient is always positive, K3 > 0, but we can have either K2 > (extensile 'pusher'), K2 = 
(symmetric quadrupolar swimmer) or K2 < (contractile 'puller'), depending on the choice of the oscillation 
amplitudes £,r^f- It is important to note, however, that the geometric restriction ^ > 2 A + S.R. + ^f implies that 
1^2 ^ K3. Consequently, the quadrupolar Ka-term effectively dominates the flow field (j27p unless one considers 
unreasonably large distances [Ij]. This is illustrated in Fig. [7] which shows the normalized flow field v := v/\v\ from 
Eq. ([27)1 and also its projection onto the unit vector r{t) for K2 > 0, K2 = 0, and K2 < 0, respectively. 

Tracer motion in the deterministic limit.- Based on Eq. (j27p it is again possible to derive an analytic estimate 
for the mean tracer trajectory. To this end, we consider initial conditions with n = (1,0,0) as depicted in Fig. [1] and 
make use of the approximation (1211) , which holds in the far-field scattering limit where Eq. (j27|) is valid. Integrating 
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FIG. 8: Deterministic tracer scattering by a ( quasi- )quadrupolar 3-sphere swimmer. In the deterministic limit Do = 0, the 
tracer trajectories converge to closed loops if the force-free 3-sphere swimmer starts sufflciently far away from the tracer, 
Xo — >■ — oo. The dotted lines in either diagram represent the analytical estimate from Eq. (|29p . (a) Results from the stroke- 
averaged dynamics for the time interval [0,55s]. The tracer particle starts at a;(0) = and the swimmer at Xq = — 150A, 
Vb = 10^. Symbols were obtained by numerically integrating Eq. (|2a|l for Do = using the averaged flow field (|27p . (b) Results 
from the corresponding time-resolved simulation for the same time interval and same initial conditions for the geometric center 
of the swimmer. The tracer oscillates due to the periodicity of the swimming stroke. Swimmer parameters: I — 10/im, A = 
0.2£,Lo — lO^s"^, A(j!) — 7r/2, which results in V ~ 12/im/s, K2 — ±0.04, K3 ~ 22. Tracer radius: ao — l/im. Simulation time 
step: At = 2.2 x lO^'^s. 



the equations of motion (|2ap for Dq 
by x{t) = {x{t),y{t),0), where 



0, we find that the trajectory of the tracer particle in the (z = 0)-plane is given 
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Note that the dipolar K2-contributions are formally equivalent to the result ([22]) for the 2-sphcre swimmer, but with 
V and K2 now determined by the geometry of the 3-sphere swimmer. Letting t — > 00 we find that asymptotically 
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which implies that the tracer particle is displaced by a finite distance due to the effective short-range interaction with 
a >2. Equations ((30|) predict that, if the swimmer starts very far away, corresponding to the limit Xq — ^ — cx), the 
tracer returns to its initial condition x[0) = 0, thus forming a closed loop in the (z = 0)-plane. This is similar to our 
previous result (j23p for the 2-sphere dipolar swimmer. However, comparing with Fig. [5] (a), the shape of the loop in 
Fig. m (a) differs significantly, because the flow fleld (P7| of the 3-sphere swimmer is dominated by the quadrupolar 
K3-contribution (for the initial conditions considered). 

The analytical estimate (|29p is based on the stroke-averaged fleld and it is interesting to compare with simulations 
for the corresponding time- resolved swimmer. The results for the latter case are depicted in Fig. \8\ (b). Due to 
the oscillations of the swimmer's legs, the tracer particle performs an oscillatory motion around its mean trajectory. 
The relative magnitude of the oscillations is bigger than in the dipolar 2-sphere swimmer case, since the effective 
quadrupolar ffow fleld of the 3-sphere swimmer decays more rapidly, resulting in a smaller absolute loop-size at 
similar initial conditions and swimmer speeds. As evident from Figs. HI the stroke-averaged result (P5)) captures the 
main features of the mean tracer motion in the time-resolved flow fleld but there is a quantitative difference of about 
20 to 30 percent. However, we observe that, asymptotically, the time-resolved tracer trajectory approximately returns 
to its initial position, if the swimmer starts sufficiently far away from the tracer particle. 

Brownian motion effects.- When Brownian motion is taken into account, corresponding to Dq > in Eq. (j2ap . 
individual tracer trajectories may differ strongly from each other. As for the 2-sphere dipolar swimmer, loop-shaped 
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FIG. 9: Brownian tracer scattering by a (quasi-) quadrupolar 3-sphere swimmer. For Do > loop-like scattering patterns can 
be recovered by averaging over many sample trajectories with identical initial conditions a3(0) — 0. (a) Results for the stroke- 
averaged dynamics obtained by numerically integrating Eq. (|2ap with Do = 0.22/im^/s. (b) Results from the corresponding 
time-resolved simulations. The filled symbols in both figures were obtained by averaging over 1000 sample trajectories, unfilled 
symbols show the results of simulations with Do — 0, and dotted lines represent the analytic estimate from Eq. (|29p . Statistical 
error bars correspond to the sample standard deviation, divided by 10 x y/i for better visibility. Prior to rescaling the length of 
the error bars agreed with the theoretically expected value \/2Dot. Compared to Fig. |8l we used the same swimmer parameters 
but the swimmer was started nearer to the tracer at Xo = —10^4 and sample trajectories were recorded for the shorter interval 
t £ [0, lis]. This explains why the loops are not fully closed. 



scattering patterns can be reconstructed by averaging over a tracer ensemble with identical initial conditions. Figures[H] 
(a) and (b) show mean tracer trajectories obtained by averaging over 1000 sample trajectories. We used the same 
swimmer parameters as in Figs. [3 but the swimmer was started closer to the tracer at Xq = —10^. The simulation 
time per run was limited to lis, as this corresponds roughly to the time scale during which an isolated alga performs a 
quasilinear motion 7]. The statistical error bars in Figs. [9] illustrate the sample standard deviation and were divided 
by 10 X -s/t for better visibility. Prior to this rescaling the length of the error bars agreed well with the value y/2Dot 
expected theoretically. In both Fig.|9](a) and (b), we notice a small, systematic drift of the mean Brownian trajectories 
compared to the corresponding deterministic trajectories (unfilled symbols) and the analytic estimate (dotted lines). 
The shift occurs in the direction of the swimmer motion and is probably caused by a noise-induced Stokes drift [20| 
and an additional, noise-induced bias due to the gradients of the flow field. 



IV. NUMERICAL METHODS 



One of our objectives was to compare the tracer dynamics in the stroke-averaged and explicitly time-dependent flow 
fields of dipolar and quadrupolar model swimmers. For this purpose, we focussed on multi-sphere swimmer models 
that can be implemented numerically in Brownian dynamics simulations [16l llTl l2ll . |22| . 

The time-resolved, 'microscopic' dynamics of the 2-sphere swimmer was computed numerically by using a spring- 
based model [17]. Within this approach, the two spheres are linked by a harmonic spring (stiffness fco) with periodically 
varying equilibrium length L{t) determined by Eq. (|16p . Compared to the implementation in Ref. jl7] , the only 
technical difference arises because our dipolar 2-sphere swimmer is characterized by an oscillating sphere radius, 
leading to time-dependent Stokes friction coefficients %{t) = 67rryAi(i), i = 1,2. If the stiffness parameter fco of the 
spring model is chosen large enough, the swimmer essentially behaves like a shape-driven swimmer (in our simulations 
we fixed fco = 0.001 kg/s^). 

The time- resolved microscopic dynamics of the linear 3-sphere swimmer can be implemented in a similar manner [l6| . 
Analogous to the 2-sphere swimmer case, two neighboring spheres (constant radius A) are connected by a time- 
dependent harmonic spring (the 'leg') with the same value of the stiffness parameter fco as above. Each spring has 
an oscillating equilibrium length given by Eq. (I24p . The velocity field acting on the tracer and the hydrodynamic 
interactions between the spheres were calculated using the Rotne-Prager-Yamakawa-Mazur tensor [23l - |27| . 

Upon adopting the same swimmer parameters for both the microscopic and the averaged model, one finds that 
the average swimmer speed measured in the microscopic simulations agrees well with the analytical estimates from 
Eqs. (ITO)) and (P5|) . respectively. We also compared the Lagrangian net fiuid flow, which was obtained from the time- 
resolved simulations by tracking a grid of non- interacting tracer particles over a stroke period [t, t + 2TT/uj]. The dipolar 
and quadrupolar flow field structures obtained by this method were found to be identical to those from Eqs. ((20)) 
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and ^7} depicted in Figs. |4] and [H 

The equations of motion for both swimmer and tracer were numerically integrated with a two-step Heun integrator 
using Nvidia's CUDA toolkit. The simulations for the time-resolved model were performed on a consumer-level 
GTX 275 GPU embedded in an Intel i860 PC running Gentoo Linux. The GPU-based parallel computation of sample 
trajectories in the case of Brownian tracer particles leads to a considerable speed-up (up to a factor of a few hundreds) 
compared with a conventional un-parallelized CPU-based implementation [28|. 

V. SUMMARY & DISCUSSION 

We have studied the hydrodynamic scattering of a tracer particle by different types of model swimmers. We showed, 
both analytically and numerically, that scattering by a force- free swimmer may lead to quasi-closed, loop-shaped tracer 
trajectories. The shapes and orientations of the loops are a signature of the properties of the stroke-averaged flow 
field. A detailed comparison with the time-resolved scattering dynamics implies that a time-averaged description 
is sufficient to capture the main features of swimmer-tracer scattering. Our analysis has focussed on two specific 
2-sphere and 3-sphere swimmer models that are analytically tractable and can be implemented numerically. However, 
the result (|29)1 is generic; it holds for the large class of low Reynolds number swimmers that generate a dipolar or 
quadrupolar flow as given in Eq. (j27p . 

We have also compared the scattering due to swimmers to that resulting when a colloid is pulled past a tracer. 
This allowed us to demonstrate the role played by the range of the interactions between the scatterer and the tracer. 
For interactions decaying as ^ l^"!"", the tracer is slowly dragged to infinity by the particle for the colloidal case 
a = 1 (within the zero Reynolds number approximation), whereas its displacement remains finite for scattering by a 
force- free swimmer with a > 2. 

We conclude the paper by addressing potential experimental realizations and by commenting on generalizations. 
The 3-sphere swimmer considered in the last part of the paper could possibly be realized in practice by adapting 
the experimental setup of Leoni et al. [8], who constructed a colloidal linear 3-sphere pump using optical tweezers, 
to mimic a self-motile, force-free swimmer. With regard to hydrodynamic scattering experiments, an advantage of 
artificially created colloidal microswimmers lies in the fact that one could tune the fluid viscosity to suppress Brownian 
motion effects in the tracer dynamics. In the case of a colloidal 3-sphere swimmer one should expect to observe loops 
that are shaped similar to those in Fig. El 

In general, Brownian motion can obscure individual tracer trajectories, but our analysis shows that averaging over 
many trajectories may allow one to reconstruct the loops even from noisy data. This suggests that such patterns could 
also be studied in biophysical experiments similar to those of Leptos et al. ^. These authors investigated how the dif- 
fusion of (an ensemble of) tracer particles in water is affected by the presence of self-motile unicellular Chlamydomonas 
reinhardtii algae. Indeed loop-shaped tracer trajectories (or, at least, segments thereof) have already been observed 
in their experiments, see inset of Fig. 1(c) in Ref. 6]. By systematically analyzing the trajectories of individual tracer 
particles in a sufficiently large, dilute suspension of algae, one could achieve the required sample size to reconstruct 
the exact structure of the loops and thus infer details about the fluid flow generated by the microorganisms. 

Lastly it is plausible to assume that an improved understanding of the underlying, elementary swimmer-tracer 
scattering processes may eventually help to model quantitatively the anomalous diffusive behavior of tracer particles 
in dilute, active suspensions (6|. The aim would be to derive an analytic expression for the time-dependent position 
probability density function of a tracer particle by averaging a representative trajectory, resulting from a superposition 
of the flow field contributions from many algae, over a set of suitably chosen random initial configurations. However, 
the averaging procedure is both mathematically and computationally demanding; we hope to tackle this problem in 
future work. 
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